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Abstract 

We show how the resummation of large logarithms can be incorporated 
into the method of effective charges. As an example, we apply this approach 
to the event shape variables thrust and heavy-jet mass in e + e~ annihilation. 
We find that, although the resummation creates problems with the behaviour 
of the effective charge in the 2-jet limit, smaller power corrections are required 
to fit the data compared to the standard approach. In addition, increasing 
the logarithmic accuracy reduces the size of the power corrections further. 
We also consider "predicting" the sub-leading logarithms in the MS scheme, 
obtaining surprisingly good results for the first few NLL terms. 
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1 Introduction 



Over the past decade, QCD has moved on from giving a qualitatively good 
description of strong interaction physics and entered an era of real quanti- 
tative tests. A wealth of precision data has appeared, particularly from the 
experiments at LEP and HERA, and it is a challenge for theorists to attain 
a similar precision in their calculations. A major stumbling block in this 
attempt is the computational complexity of QCD perturbation theory. In 
practice this means that we are limited to using the first 2 to 3 terms of 
the perturbation series in the strong coupling as- The relatively large value 
of as compared, say, to oqed, means that this unavoidable truncation of 
the perturbation series seriously limits the accuracy of perturbative QCD 
(pQCD) calculations. This problem is further complicated by the spurious 
dependence of fixed-order results on the renormalization scheme (RS) and 
scale (/i), which together can be termed the renormalization prescription 
(RP) [1 j. This dependence, which does not afflict the exact results, requires 
us to choose some suitable RP to make pQCD predictions. There are in prin- 
ciple infinitely many ways to do this; by far the most common in practice is 
to choose the MS RS and set /i equal to some physical scale of the process. 
For brevity, let us refer to this approach as MSPS (for MS scheme with a 
physical scale). Despite its popularity, this choice of RP is not supported by 
any theoretical argument; the motivation for it is purely empirical, namely 
the evident success of standard pQCD phenomenology, where this approach 
is almost invariably taken. However, there also exist several theoretical pro- 
posals for choosing so-called "optimized" RPs, tuned to the observable at 
hand using information that is already present in the perturbative calcula- 
tions. This is in contrast to the idea of selecting a "physical" scale which 
requires an additional judgement to be made for each observable outside the 
well-defined framework of pQCD. Of course, the choice to use a particular 
"optimized" approach is also an additional judgement, but the point is that 
these approaches are supposed to pick the RP for each observable automati- 
cally once the optimization method itself is specified. 

One such optimization method, the method of effective charges (ECH) 
was recently applied in Ref. [3j to the description of event shape means. 
The authors found that with this choice of RP, next-to-leading order (NLO) 
pQCD could describe the data very well even without the introduction of 
power-suppressed non-perturbative effects. This surprising result suggests 
that it might be interesting to apply ECH also to the distributions of these 
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event shapes. Such analyses have been performed |U |S] , but unfortunately 
in this case ECH suffers a restricted range of validity due to kinematical 
end-points and logarithmic enhancements. This latter problem affects also 
the MSPS treatment of event shape distributions, where it can be handled 
by a resummation of large logs. A way to implement this within ECH was 
presented in Ref.jU], but a systematic recipe for carrying out the ECH analysis 
of a general observable known to, say, next-to-leading logarithmic accuracy, 
is not available in the literature. In this paper we present such a recipe. We 
then go on to test it by applying it to the problem of describing event shape 
distributions in the 2-jet limit. 

The plan of this paper is as follows: we briefly review ECH in Section|21 In 
Section|H]we outline a general procedure for obtaining an all-orders resummed 
ECH result starting from an all-orders resummed MSPS perturbation series. 
Section |U gives our application of this to the distributions of the event shape 
variables thrust and heavy-jet mass in e + e~ annihilation. Section El contains 
our conclusions. 



This section briefly summarises the method of effective charges; for more 
detailed developments see e.g Refs. [3 E] . Consider a quantity normalized so 
that its perturbation series takes the form 



a(/x, RS) + nifx/Q, RS)a 2 0u, RS) + r 2 (/i/Q, RS)a 3 (//, RS) + • • • , (1) 



with a = ag/ft- Such a quantity could either be the coupling defined in some 
RS with Q as the renormalization scale or a suitably normalized physical 
observable H(Q) depending on a single energy scale Q. Indeed, these possi- 
bilities are not mutually exclusive. Given some observable of the form Eq. 
(JTJ), one can define an RP such that all the r n vanish; in this RP, 1Z = a, so 
the observable is equal to the coupling. Such an observable/coupling is called 
an effective charge j2]. A typical example is the R-ratio in e + e~ annihilation, 
with which one can associate an effective charge TZ(Q) via 



\ / / 

Whatever the quantity represented in Eq. (0), it will be independent of 
the RP in which the expansion is performed. Moreover, when the expansion 
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is truncated at some order, the variation in this partial sum due to a change 
of RP is always one order higher in a. This implies specific relations between 
the r n and the RP 0. 

The dependence of the coupling a defined in some RS on p is described 
by the beta-function 

^dln^ = ^ = ~ ba2 ^ + ca + c 2° 2 + c 3« 3 + •••)> (3) 

where b and c are independent of the RS, and the q, i > 2 can be taken to 
label the RS, along with the scale parameter A which we introduce below jSJ. 
Restricting a to the subset of couplings which are also effective charges, the 
same equation is conventionally written as 

= P {u) = -bn\i + cK + P2 n 2 + P3 n* + •••), (4) 

where b and c are now independent of the choice of effective charge, whereas 
the p n depend on this choice. 

Because p describes a functional relation between two physical quanti- 
ties, namely 7Z and its energy derivative, it is automatically independent 
of any choice of RP we make in calculating it. Clearly, the p n must also 
be similarly RP-independent. This means that truncating p at some fixed 
order in perturbation theory, and calculating an approximation for 7Z by in- 
tegrating the corresponding approximate effective charge beta-function will 
give RP-independent results. Another way of looking at this is to say that 
the method of effective charges involves a specific choice of RP (i.e. the RP 
where TZ = a), so that the energy evolution of the observable is identical to 
the beta-function evolution of the coupling. 

Now, to obtain a perturbative approximation to p, consider setting p = Q 
in Eq. (0) (which is allowable as the full sum is independent of p) , and then 
differentiating with respect to InQ. This gives 



- dn d1Zda(Q) 

P\/'C) 



= -~-P{a) 

a=a(K) OCl 



(5) 

a=a(TZ) 



dhiQ da din Q 

Expanding this order-by-order in TZ and comparing to Eq. (J3J) gives expres- 
sions for the p n as multinomials in the r n and the c n ; for example, the first 
two are 

P2 = c 2 + r 2 - T\C - r\ (6) 
P3 = c 3 + 2r 3 - 4rxr 2 - 2r x p 2 - r\c + 2r\. 
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Suppose we have performed a NLO calculation of 1Z (in some arbitrary 
RP). This provides us with r 1; but not r 2 or any higher-order coefficients. 
So as this stage, our best approximation to p(TZ) is simply the universal 
form p(7Z) = —bTZ 2 (l + cR). This can therefore be termed the "NLO" 
approximation to p. Given the and c, for i = l..n, we can calculate the pi 
to the same order, obtaining a N"LO approximation. 

Having constructed an approximation for p it just remains to relate this 
to the effective charge itself and hence to the corresponding observable. This 
can be accomplished by integrating Eq. (JU) using asymptotic freedom as 
a boundary condition. Note however that Eq. (£Q) only fixes 1Z up to an 
arbitrary rescaling of Q. This is necessary, as the NLO p is the same for 
all effective charges, and they certainly should not all have the same NLO 
predictions. Similarly, Eq. (j3J) only fixes a up to a rescaling of p. 

The missing information needed to fix the scale (corresponding to the 
single free parameter of massless QCD) can be conveniently provided by 
a parameter A associated to each RS and each effective charge which will 
arise as a constant of integration. It might seem contradictory to describe 
A as simultaneously differentiating between effective charges and acting as 
the free parameter of QCD. The resolution of this apparent paradox is the 
so-called Celmaster-Gonsalves relation which relates A parameters for 
different effective charges (or indeed different RS's) exactly using only their 
NLO coefficients computed in some reference RP 



(It does not matter what RP is used for the calculation of these r% coefficients 
as a change of RP just adds a constant to r 1 which clearly drops out of 
this equation.) We can therefore take any of these A's as the fundamental 
(dimensional transmutation) parameter of QCD and obtain any other A from 
it given only NLO information. Conventionally we choose A^jg. However, 
there is a subtlety here related to the way that Aj^g is traditionally defined, 
so that to avoid confusion we refer to the A^jg we are discussing here as A^jg 
jH]. This is related to the more conventional A^jg [TO, via 



Ane~ ri/b = A w e 



r'Jb 



(7) 
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After integrating Eq. (jlj) we arrive at 

Q f°° dx [ n{Q) , 
bin-— = / -r- -+ dx 

A n Jll{Q) x (1 + cx ) Jo 



-1 1 

+ 



X 2 (l + CX + P2^ 2 + • • •) X 2 (l + CX) 

(9)' 

The first integral is independent of the particular effective charge under con- 
sideration, whereas the second is not. Exponentiating gives 

^ = r(K(Q))Q(W)) (io) 

where T is universal and comes from the first integral in Eq. Q 



F(K)=e- 1 "*\l + —) , (11) 
and Q comes from the second integral and depends on the effective charge, 



Q(TZ) = exp 



dx—— + 



o p(x) bx 2 (l + cx) 



(12) 



Approximating p by its NLO form gives {?nlo — 1- Converting from An to 
Aj^g using Eq. (JZj) and Eq. (jHJ gives 

A m = Q^(n(Q))g(nQ))e- ri{m " =Q)/b (jj ' ■ (13) 

This equation allows us to extract values of A^jg directly from the observed 
values of 7Z, or to make predictions for 1Z (by solving the implicit equation 
e.g. iteratively). 

One way to compare this with the more standard approach of truncating 
the series for TZ in some fixed RP is to write an "effective" effective charge 
beta-function describing the energy evolution of 1Z implied by this standard 
approach. Let us call this function p. It can be computed from Eq. (J3J) using 
the truncated relation between 1Z and a. It always agrees with the true p 
up to the order (NLO, NNLO,...) to which 1Z has been calculated, but it 
also includes terms of all higher orders in TZ which depend on the RP we 
chose to perform our truncation in. If the RP is such that r\ is large, these 
higher order terms are also large, so the ECH predictions will differ radically 
from the standard prediction in that RP. For many e + e~ jet observables T\ 
is indeed large in the MS PS RP. 
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3 Resummation of Logarithms in the ECH 
Beta- Function 



It is commonly stated that the method of effective charges is inapplicable 
to exclusive quantities such as event shape distributions. The idea is that 
the dependence of the physical quantity on multiple scales invalidates the 
derivation of the ECH beta-function as presented here in Section EJ How- 
ever, as pointed out in [7j, this is not really the case. Given an observable 
TZ = TZ(Qi,Q 2 , ■■■,Qn) depending on n scales, one can simply re-express it 
as TZ = U{Qx,Q 2 /Qi, -,Q n /Qi) = T^x 2 ,...,x n (Qi)- Here the = Qi/Qi 
are dimensionless quantities that can be thought of as labelling the effective 
charge which is now a function of one single dimensionful scale Q\. We can 
then write an effective charge beta-function for this TZ describing the energy 
evolution of our observable for fixed values of the ratios Xi. Still, this for- 
mal manipulation cannot tell us whether the p function we arrive at in this 
way will be well approximated by its NLO terms, which is what we require 
for most fixed-order phenomenological applications, given the current state 
of the art in perturbative QCD calculations. One reason in particular why 
this might not be the case is if some of the Xi become large - typically this 
leads to powers of large logs = log(xj) enhancing the coefficients r n in 
the perturbation series for TZ, and hence also the p n in the corresponding 
p{TZ) function. A common situation is that more logs appear as the order of 
perturbation theory is increased, so that for a sufficiently large L the terms 
all become of similar magnitude. This invalidates both the NLO (universal) 
approximation for p(TZ) and NLO MSPS. However, in the latter case, there 
is a well-known way out. If the leading powers of the logs can be identified 
as having some simple dynamical origin, we may be able to calculate them 
to all-orders and then effect a resummation, extending the validity of our 
perturbative results into the large L region. This suggests that essentially 
the same trick might work for the p function. In this section we describe 
how to accomplish this, expanding on ideas in Ref. jHj; an example of the 
phenomenological application of these ideas to event shapes is presented in 
the next section. 

Our approach will be to start with some resummed result for the observ- 
able of interest calculated by conventional means. As an example, consider 
an observable of the form 

0(L) = LA LL (aL) + A NLL (aL) + aA NNLL (aL) + • • • , (14) 
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where a is as usual the couplant and L is the large logarithm which this 
expression resums. The subscripts "LL", "NLL" and "NNLL" stand for 
leading log, next-to-leading log and next-to-next-to-leading log respectively. 
We can relate this observable to an effective charge 

0(L) = r (L)TZ(L) = r (L)(a + ri {L)a 2 + r 2 (L)a 3 + ■■■). (15) 

Here tq is the leading order coefficient, whose large L behaviour is tq ~ L 2 . 
The r n can now be expanded in powers of the large log L. In this case their 
leading behaviour is r n ~ L n and we can write 

r n = r^L n + rl^L n - x + ■ ■ • , (16) 

so the structure of the p n as illustrated in Eq.© implies that 

Pn = p ^L n + p^ L L n - l + ..- . (17) 

Because L is a logarithm of a physical quantity, this expansion of the p n is RP- 
independent We can thus define resummed, RP-independent approximations 
to p(1Z) such as 

oo 

pll(^) = -bn 2 (i + c n + Y,Pn L L n n n ) (is) 

n=2 

oo 

Pnll(^) = -b^il + cn + Y.iP^^ + P^L^TZ") . (19) 

n=2 

and so on. These can be calculated order-by-order using the relations between 
the p n and r n , the first few of which are shown in Eq.©. Alternatively, we 
can apply a numerical procedure to extract our desired p function from 1Z 
calculated to similar logarithmic accuracy. This is particularly simple for 
Pll, as can be seen by considering the p corresponding to 7^ll with one- loop 
beta-function 

PL L (*)=/3(a)^ = -6a 2 ^, (20) 
da da 

with a chosen such that 7^ll(g0 = x. The perturbative coefficients of this p 
function can be obtained from the expressions for the p n , of which the first 
two are shown in Eq.©, using c = 0, q = and r n = r^ L L n . But then it is 
easy to see that the coefficients we obtain are proportional to L n and moreover 
identical to the coefficients of pll (because adding the sub-leading terms in 
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(3(a) and 1Z only affects the p n at next-to- leading logarithmic accuracy). In 
other words /5 L l, defined as above, is equal to pll except for the clZ term 
which can easily be added (in some sense this term is NLL as it is 0(L~ 1 ) in 
the large L, fixed 7ZL limit, but we include it in pll as it is obviously present 
in the full expression, and this avoids having to modify Eq. (fT2|) ). Thus 
Pll(TZ) can be calculated to arbitrary accuracy for a given 1Z by numerically 
inverting TZ-^^a) to obtain the corresponding a, then substituting this a into 
Eq.®. 

Going beyond leading log accuracy things become slightly more compli- 
cated, because the p functions pick up terms at lower logarithmic accuracy 
that do not appear in our resummed approximations. For example, at NLL 
we will have r\ = r\ L L + r^ LL , so that p2 will contain not only L 2 and L 1 
terms, but also L° terms. These are not included in our definition of pnll, 
and indeed they must not be, as they are affected by the addition of the re- 
maining missing terms in 7Z and hence are RP-dependent. However, it is still 
the case that the NLL terms in Pnll are unchanged by adding sub-leading 
terms (in (3(a) and 7Z), and are therefore identical to the corresponding terms 
in Pnll (assuming that the p functions are defined with beta-functions hav- 
ing sufficiently many terms to make this true, e.g. (5(a) = —ba 2 (l + ca) for 
the NLL case). So, truncating p by numerically taking limits (L — > oo with 
LIZ fixed) allows us to extract the LL and NLL terms. The generalization 
to higher logarithmic accuracy is straightforward. 

Some physical quantities might have more divergent logarithmic behaviour, 

eg. 

O(L) = A LL (aL 2 ) + L~ l A NLL (aL 2 ) + L~ 2 A NNLL (aL 2 ) + ■■■ . (21) 

In this case, r n ~ L 2n but the preceding argument goes through essentially 
unchanged, except that p n ~ L 2n as well. 

Having obtained a resummed p function, we can proceed to extract A 
from the observable making use of Eq.Q with 1Z = O(L)/r (L). In doing 
this, we may have available exact values for r and r x from a fixed-order calcu- 
lation, which we can use in place of their approximations from the resummed 
results. This allows us to combine resummed and fixed-order information in 
an essentially unique way (once we have fixed the definition of our effective 
charge), avoiding the so-called "matching ambiguity" associated with doing 
this in MSPS. In particular the full exact NLO coefficient r\ in a given RS is 
reproduced if 1Z which solves Eq. © is expanded in the coupling a for that 
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scheme, thanks to the RS invariant A-^ which appears on the lefthand side 
of the equation. 

This provides us with a complete method for making ECH predictions 
including resummations of large logarithms. In the next section we test this 
approach by comparing the distributions of thrust and heavy-jet mass in 
e + e~ annihilation to NLL ECH predictions. 

4 ECH for Event Shapes at Next-to-leading 
Logarithmic Accuracy 

Event shape variables provide some of the most interesting and useful ways 
to confront QCD calculations with experiment (for a recent review see [11J). 
Their infrared and collinear safety guarantees that they can be calculated 
in QCD perturbation theory, but fixed order calculations describe their dis- 
tributions rather poorly. This situation can be improved by the recognition 
that for a shape variable y, vanishing in the 2-jet limit, large logarithms 
L = log(l/y) appear at each order of perturbation theory and must be re- 
summed. However, a full understanding of the observed distributions re- 
quires the introduction of large non-perturbative effects, power corrections 
oc e~ h l a ~ A/Q where Q is some hard scale (e.g. the e + e~ centre of mass 
energy). Although the existence of such effects can be motivated by consid- 
ering simple models of hadronization o r through a renormalon analysis 
[T5j . their magnitudes are not at present calculable in a truly systematic 
way from the QCD Lagrangian. Therefore to fit the data we require the 
introduction of either a phenomenological hadronization model or additional 
non-perturbative parameters. Although this hampers attempts to extract 
reliable measurements of as from the data, it also provides a good oppor- 
tunity to study the IR behaviour of QCD experimentally. For example, in 
Ref. [H] it was proposed to relate the magnitude of the 1/Q power correction 
to event shape means to the average value of a hypothetical infrared-finite 
coupling. This approach leads to predictions for 1/Q power corrections to 
all event shape means in terms of a single additional parameter, a Q (/ii), 
the zeroth moment (i.e. mean) of the coupling at scales < // < [ii (typ- 
ically hi ~ 2GeV). In Ref. [To] it was shown how this approach could be 
extended to apply to event shape distrubutions. Since then many experi- 
mental studies have appeared, fitting event shape means and distributions 
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simultaneuosly for a^(Mz) and a (fij). Generally an approximate (up to 
corrections ~ 25%) universality of the a values is observed, supporting the 
hypothesis that power corrections can be related to a universal coupling in 
this way However, most of these fits use the MS scheme with the "phys- 
ical" scale choice \i = Q. An exception is work within the Dressed Gluon 
Exponentiation (DGE) framework (see Refs. jTol IT7| ITS] for the application 
to average thrust, and the thrust and heavy-jet mass distributions). For the 
event shape means another possible approach is to work within the ECH 
framework. This was first carried out for 1-thrust in Ref . [TTi] . and some- 
what reduced power corrections were found compared to the physical scale 
approach. This suggests the possibility that the power suppressed effects are 
partly compensating for missing higher-order perturbative terms associated 
with the running of the coupling. Indeed, similar conclusions were reached in 
the DGE approach, although there the subset of terms resummed differs sig- 
nificantly from that resummed by the change of scale implicit in NLO ECH: 
the former is factorially divergent, the latter actually converges. Recently an 
analysis similar to that of Ref.|19j. though more extensive, was performed 
by the DELPHI collaboration 3J, taking into account effects arising from 
the finite bottom quark mass via Monte Carlo simulations. They found re- 
markably small power corrections within the ECH approach, which for many 
observables were consistent with zero. Indeed, the ECH predictions with no 
power corrections whatsoever gave a better description of the data than the 
model of Refs. O HH] with a universal a - In light of these surprising results, 
it is interesting to consider applying the ECH method to the event shape 
distributions. 

In fact, event shape distributions have previously been studied within 
the ECH framework Ref. [3] in particular studied how the fit of the 

ECH results to data varied in quality for different regions of phase space. 
To do this an effective charge was constructed separately for each bin of 
the data, and NLO QCD calculations were used to extract A^g at centre 
of mass energy Q = M Z - The consistency of these Aj^g values between 
different data bins could then be examined. Non-perturbative effects were 
taken into account by using Monte Carlo simulations to correct the data back 
to "parton-level" distributions. This generally improved the consistency of 
the Ajyfg measurement, but with this approach it is hard to see whether the 
ECH distributions prefer smaller hadronization corrections than the MSPS 
ones. Moreover, even after these corrections were applied there were still two 
kinematical regions where the effective charge ceased to be a good description 
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of the data, leading to instability in the A^g values: the 2-jet limit where 
large logs enhance the higher-order perturbative coefficients, and the region 
(which exists for many observables) where the LO result vanishes, causing 
7*1 — » oo. The latter problem is hard to address within the effective charge 
approach, but the former problem seems to call for a resummation of the 
effective charge beta-function, as described in Section 01 

In this section we study the distributions of 1-thrust (r = 1 — T) and 
heavy-jet mass (ph)- We first show the effect of replacing the hadronization 
corrections of j3] with an analytical power correction ansatz. For simplicity, 
we use a shift in the distribution by an amount C\/Q. This form can be 
motivated by considering simple models of hadronization or through a renor- 
malon analysis ^2] and has been found successful phenomenologically (see 
for example p.). Although better fits are often obtained using the model 
of [TU [T5] , because we are using a different perturbative approximation to 
standard NLO QCD, the subtractions needed to remove double counting will 
not in general be so simple - in particular, it is not clear what scheme should 
be used for performing the subtraction. 

Next we consider placing the effective charge into the exponent of the inte- 
grated thrust distribution. Even using a NLO approximation for the effective 
charge, this has the effect of resumming a series of logs in the distribution 
itself (in particular the "double logs" are included). 

Finally we present results showing the effect of using the resummed p 
functions described here in Section |H1 First we investigate the extent to 
which higher order MSPS logs are already included in the lower order ECH 
predictions ("RG-predictability"). Then we actually perform fits using the 
resummed ECH predictions. 

For comparison, at all stages we also give results of fits to the same data 
using MSPS (at NLO, LL and NLL accuracy). As is customary, we use the 
variation of p such that Q/2 < p < 2Q to estimate a "theoretical error". 

The general question of separating perturbative and non-perturbative ef- 
fects also deserves comment. Because perturbation theory diverges, it is not 
straightforward to define its sum; however, without doing this the magnitude 
of the "non-perturbative" effects is ambiguous. Therefore, it is preferable to 
combine a fit for power corrections with a renormalon resummation, as in 
[To! IT7| ITH] . We have not done so in this analysis, but as we are comparing 
ECH to MSPS which differ only by a convergent set of higher order terms, we 
believe that our conclusions regarding the relative size of power corrections 
stand. It would of course be interesting to investigate the effects of a renor- 
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malon resummation on our ECH results (ECH renormalon resummations 
have already been carried out for some single-scale observables in |2Uj). 

Our data is taken over a wide range of centre-of-mass energies Q = 35 — 
189GeV fRefs.|5].|21|-|31j). Lacking information on the correlation between 
data points we have simply combined statistical and systematic errors in 
quadrature and performed a min-% 2 fit, allowing x 2 t° vary by 4 from its 
minimum to estimate a 2a error. This over-simplistic treatment means that 
our errors cannot be considered reliable, however the central values of A-jjg 
do give an impression of the effect of including the power corrections and 
logarithmic resummation into the ECH framework. 

Let us now consider the effects of analytical power corrections on the 
results of The procedure used in |3J was to write an effective charge to 
represent the value of the event shape distribution integrated over each bin of 
the data. First, the Monte Carlo program EERAD [32] was used to compute 
the NLO perturbative coefficients 3 for each bin 



These were then used to write an effective charge, from which a value for 
Aj^g could be extracted by feeding the data into Eq. (fl3~j). Here, to introduce 
a fit for C\ we drop this "direct extraction" approach and instead perform a 
minimum y 2 fit f° r A^jg and C\. For this to work, we need to exclude the 
regions where the EC approach cannot fit the data. For comparison with 

we choose the same ranges selected there (based on the flatness of ri), 
except that the lower end of the range is made proportional to 1/Q when 
looking at data away from Q = Mz- The reason for this is that sub- leading 
non-perturbative effects are expected to become important for y ~ A/Q 
[TBI IHH ES]- As we are shifting the predictions before comparing to data 
we require the NLO coefficients evaluated for arbitrary bin edges. We have 
approximated these using a set of order 6 polynomial interpolations from the 
output of EVENT2. We have checked by halving the Monte Carlo bin size to 
0.005 that this induces no sizeable error (using the doubled sampling changes 
the best fit values here by less than 2%). 

The results for 1-thrust and heavy-jet mass are presented in Fig. 1 and 
Fig. 2 respectively. 

3 For our analysis we actually used EVENT2 J2] and we have checked that both pro- 
grams give consistent coefficients. 
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Figure 1: 1-thrust: Fits for A^g and C\ within the framework of [3] (solid 
ellipse) and standard NLO QCD peturbation theory (dashed ellipse). In the 
latter case the scale is chosen to be fi = Q, and the effect on the central 
value of a change of renormalization scale by a factor of 2 is indicated by the 
arrows. 2a error contours are shown (from allowing \ 2 to vary within 4 of 
its minimum). The fit range is 1 — T = 0.055Mz/Q — 0.23. 
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Figure 2: As Fig. 1 but for heavy-jet mass. The fit range is 
0.035M Z /Q - 0.2. 
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In the case of thrust it appears that the ECH results prefer larger power 
corrections (and significantly smaller Aj^g values). For heavy-jet mass, the 
situation is similar, although the differences are not quite as extreme. How- 
ever, in both cases we find comparable A^jg values to those found in jl] us- 
ing hadronization corrected data rather than an analytical power correction 
ansatz. 4 

A crucial property that an event shape must possess in order for a re- 
summation of logarithms to be feasible with present techniques is so-called 
exponentiation. To illustrate this property, consider the typical form of an 
event shape distribution as a double expansion in a and L = \og(l/y): 

= A LL (aL 2 ) + L- l A NLL (aL 2 ) + L~ 2 A NNLL (aL 2 ) + ■■■ . (23) 

a dy 

The A functions have a perturbative expansion A(x) = Aqx + A\x 2 + . . . and 
for r and pn are known up to NNLL accuracy. If the event shape exponenti- 
ates, then we can also write 




= C(a7i) exp(Lgi(cnrL) + g^icmV) + ags(anL) + 



ag 3 (airL) + ---) + D(air,y) . (24) 

For t and ph, gi and g>2 are known |SZ|- When working with this form of 
the distribution it is conventional to refer to g\ as containing the leading 
logarithms and #2 as containing the next-to-leading logarithms. C = l + 0(a) 
is independent of y, and D contains terms that vanish as y —>■ 0. These 
can be calculated to NLO by comparison with fixed-order results. However, 
there is no unique way of including this fixed-order information into R y (y') 
(this is the so-called matching ambiguity). For example, it is also legitimate 
to include the C,D terms into the exponent (termed "log R matching"), as 
the difference is of order a 3 . 

In Eq. 12 Hj) . there are terms at 0(a n ) with up to 2n factors of L mul- 
tiplying them. In contrast, in the exponent of Eq. (|24jl at 0(a n ) we find 

4 However, the values of A^jg quoted in 0] are actually wrongly normalized for two 
reasons. Firstly the factor of (2c/b)^ ~ 0.85 was omitted, so the results are really 
values for A-^. Secondly the results of EERAD were normalized to the Born cross- 
section (To, whereas the data are normalized to the total cross-section a, and this was not 
taken into account. Mutiplying the EERAD perturbation series by a correction factor 

<7o/cr = 1 — a/ir -\ decreases r*i by exactly 1, increasing the extracted A^jg- values by 

e 1 / 6 . So the total correction factor to apply to the results of g] is (2c/b) < - c / b ">e 1 / b ~ 1.11. 
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no more than n + 1 factors of L. The L m a n terms with n + 1 < m < 2n 
are generated by the exponentiation. For example, including just the leading 
order, leading log term ~ L 2 a in the exponent of Eq. ()24j) leads, after the 
exponent is expanded out, to the entire set of double- logs ~ L 2n a n in Eq. 
(|2*3*|) . Ideally we would like to use this exponentiation property in our ECH 
approximation. So, let us consider the effect of defining 

R y (y') = exp(r (3/'W)) ■ ( 25 ) 

Here all the physics is encoded into a single effective charge, which is ex- 
ponentiated in its entirety. This is similar to log R matching in that if we 
re-expand r 7Z in terms of a and L the C and D functions will clearly appear 
in the exponent. However, in this approach there is no matching ambiguity 
because once we have picked the effective charge the inclusion of C and D 
is automatically determined. The function ro for thrust can be found by 
integrating the analytically known leading order 1-thrust distribution 

1 da 
a dr 



C F a (3 - 9 r - 3 r 2 + 9 r 3 + (-4 + 6 r - 6 r 2 ) log(± - 2)) 



lo 2r(r-l) 

(26) 

with the boundary condition that R T vanishes to LO for r > 1/3. The result 
is 

ro(r) = C ?(-l + Y + 3T + ^T + (l- 3r ) 0°g(l-2T)-log(r)) 
-(log(l - t) - log(r)) 2 - 2Li 2 (rhpj) ■ W 

1-Thrust and heavy-jet mass agree at LO, so the same result holds for heavy- 
jet mass. A given prediction for TZ(y) now allows us to calculate a correspond- 
ing R(y). Then, by taking the difference in R(y) across the bins in each 
experimental data set, a comparison to data can be carried out, including a 
1/Q power correction by using Rpciy) = Rpt(v — Ci/Q)- hi the remainder 
of this paper we will consider predictions for the distributions of 1-thrust 
and heavy-jet mass arising from subsituting various ECH approximations 
into Eq. 

The simplest possibility is to use a standard NLO ECH TZ(y). This only 
requires knowledge of ri, which can be easily obtained from the results of 
Monte Carlo calculations of the distributions to NLO. This NLO ECH, re- 
expanded in a(Q) and L in the MS scheme, includes terms ~ L m a n for all n 
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and all m < n + 1. Some of these terms can be compared with their exactly 
known LL and NLL counterparts allowing us to determine to what extent the 
LL and NLL terms are "RG-predictable" . In this context, RG-predict ability 
refers to the extent to which the higher order coefficients in the perturbation 
series for 1Z are already present in some lower order ECH result. Overall, the 
logs are not very well predicted by the NLO ECH results except for rather 
small n (see Figs. 3 and 4 for NLL examples). However, the exponentiation 
of the effective charge will produce further towers of logs in the distributions 
themselves. We might therefore expect these NLO ECH results to have better 
behaviour in the 2-jet region than the results of Ref. 0J. This is indeed the 
case, as can be seen in Fig. 5. It is interesting to note that exponentiation 
of a NLO MSPS series in place of 1Z produces a distribution (the dashed 
curve in Fig. 5) with a badly misplaced peak (this remains the case for any 
reasonable value of Aj^g). Therefore, the good qualitative description of the 
peak is only obtained at NLO with the use of both exponentiation and ECH 
(until we introduce non-perturbative effects). 5 

We can now consider performing a resummation of logs in the p function 
as described in Section EH First, constructing these order-by-order in TZ 
allows us to again address the question of the "RG-predictability" of the 
MSPS logs. For example, one can ask how much of the NL log at 0(a n ) 
in MSPS is included when we use the LL ECH result (of course, if we use 
the NLL ECH the full NL MSPS log is included by construction). To find 
out, we can re-expand the ECH results in terms of Oms(Q) as we did for 
the NLO ECH. The resulting coefficients are shown in Fig. 3 (for 1-thrust) 
and Fig. 4 (for heavy-jet mass) as fractions of the exact coefficients. The 
LO and NLO (n = 1,2) coefficients agree exactly as and ri have been 
used to NLL accuracy in all the predictions. There is a clear improvement 
in the prediction of the NL logs as we move from NLO ECH to LL ECH 
as one might expect. The extent to which the NL logs really are included 

5 At NLO, ECH is equivalent to a scale choice p — Qe~ ri / b . In the case of both thrust 
and heavy-jet mass r\ = bL/2 + const +• • •, so using a NLO ECH is equivalent to choosing 
M = Q*Jyf{y) where f(y) goes to a constant as y — > 0. This is interesting as a "physical 
scale" argument where one takes, for example, the heavy-jet mass rrih = y/pKQ as the 
scale would give p — Q^/pK- So these two scale-setting methods have the same leading 
Ph dependence, and only differ by the factor /. This factor is important, however, as its 
overall normalization is RS-dependent and this ensures that p, is chosen in such a way that 
we obtain the same ECH answer whatever RS we choose. 
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predicted/ exact 

* • 



Figure 3: Prediction of NL MSPS logs in the exponent for 1-thrust based 
on re-expanding lower order ECH results in a^jg. The ratio of the predicted 
NL cofficient at 0(a n ) with the corresponding exact coefficient is shown. 
Triangles are for NLO ECH, diamonds for LL ECH. As a consistency check, 
squares show the results for NLL ECH where the NL log must appear exactly. 
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Figure 4: As Fig. 3 but for heavy-jet mass. 
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Figure 5: Comparison of the 1-thrust distributions calculated using various 
NLO approximations in the 2-jet region. The solid curve arises from the 
exponentiated ECH of Eq. (|23j). The dashed curve is obtained by expanding 
this effective charge in the MS scheme with fi = Q. The dotted curve is a 
prediction in the approach of Ref. [3]. Throughout we have taken Q = Mz 
and Ajyfg = 250MeV. For comparison, DELPHI data at Q = Mz are shown. 
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in LL ECH is encouraging, as it suggests that NLL ECH might do a good 
job of including some higher order corrections that are omitted in the MSPS 
approach. 

The method described in Section|3]can now be used to produce numerical 
approximations to pll and pnll (these are not truncated at any order in 
7Z). All of our calculations were carried out using the computer algebra 
system Mathematica jSHj; allowing the use of arbitrary precision arithmetic 
in taking the L — > oo limit. These pll and pnll functions can be used 
to make predictions for 7Z, by inserting them in Eq. Q and numerically 
solving the transcendental equation. To ensure exact cancellation between 
the singularities in Eq.(|12|). for 1Z < 0.005 we use the exact series expansions 
of p{JZ) up to to order 7Z 4 (the difference is totally negligible for all values 
of L we consider) This defines what we call our "LL ECH" and "NLL ECH" 
predictions. 

To perform fits with these exponentiated effective charges we again need 
to select a fit range. After exponentiation, the problem as r — > remains, 
and in fact for thrust worsens; unfortunately this means we need to restrict 
the fits to 1 — T < 0.18, ph < 0.24 to obtain good fits in the ECH approach. 
Irrespective of the inclusion of logs, the onset of non-perturbative effects more 
complicated than a simple 1/Q shift means that we still need to impose a 
lower cut. These higher-order non-perturbative effects are expected to be of 
order A/ (Qy) so our cut should be placed at y ~ A/Q with A some infrared 
scale. One might expect the inclusion of the extra logs into p to improve the 
fit of the ECH prediction to data in the 2-jet region. Unfortunately, it turns 
out that including these logs actually worsens the behaviour of the ECH 
results. In this region the growth of r\ causes 1Z to become larger (because 
A-fc approaches Q), and this is accelerated by the addition of logs into p. In 
fact, 7Z eventually grows large enough that we encounter a branch cut in p 
which appears due to the branch cut in gi |S7j . Clearly this behaviour is 
unphysical, and must be avoided in our fits to data. Presumably some other 
higher order corrections intervene to produce an ECH prediction which is 
well-behaved in the 2-jet limit. In any case, in light of the RG-predictability 
of the sub-leading MSPS logs, it is still possible that the p resummations 
improve the quality of the ECH predictions in the intermediate region, so it 
is worth trying to fit the data with a lower cut in place (and studying the 
sensitivity of the best fit parameters to the choice of this cut). Good fits are 
obtained over the whole energy range using p h , 1 — T > 0.05M Z /Q- 

Any data bins not lying within the range 0.05Mz/Q < p^ < 0.24, 
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0.05M Z /Q < 1 — T < 0.18 have been left out of the fit; a summary of 
the data we actually used is given in Tables 3 and 4. We have also removed 
the JADE data at 35 and 44 GeV from the heavy-jet mass fits, since its 
inclusion dramatically worsens the fit quality for both the MSPS and ECH 
predictions. The result of fitting for A^jg and C 1 are shown in Figs. 6 and 7 
for all three approximations. For comparison, MSPS results are also shown. 
The fit range for these could in principle be extended as they do not suffer 
from the tq — > problem that afflicts the ECH, however, in order to facilitate 
a direct comparison between the two approaches we have used the same fit 
range for both. The most notable feature of the results is the stability of 
the Ajfg values found within the ECH framework as we move from NLO to 
LL and then to NLL accuracy, while the fit quality hardly changes. It shows 
that the effect of the leading and next-to-leading logarithms can be mimicked 
by an increase in C\ (which itself is not large). This improved stability with 
respect to the order of the approximation might be a consequence of the 
RG-predictability of the MSPS logs discussed above, because, for example, a 
lot of the logs that only turn up at NLL order in the MSPS predictions are 
included already at LL order in ECH. It must be noted however that despite 
their stability, these Ajjg values are still smaller than the world average. 
Some examples of the actual NLL ECH distributions are shown in Fig. 8. 

To investigate the sensitivity of these results to our choice of fit range we 
have redone the fits for a "low" range and a "high" range. The low range was 
determined by decreasing the upper cut until half the bins were excluded, 
and the high range was determined by increasing the lower cut similarly. 
The effects of these changes on the central values of and C\ are shown 
in Table 1 (for 1-thrust) and Table 2 (for heavy-jet mass). The ECH fit 
values for heavy-jet mass appear more stable than the MSPS ones; and the 
stability increases as the accuracy of the predictions are increased. This is 
also true for the thrust, but to a lesser extent. A particular exception is 
that the application of NLL ECH to the "high" fit range gives a significantly 
different power correction compared to the "low" fit range. This is probably 
responsible for the relatively large \ 2 f° r the "normal" fit. The reason for the 
change in C\ may be the r$ — > problem being exacerbated by the increase 
in size of the effective charge as more logs are added into its beta-function. 
Because the r — > problem represents a breakdown of our approximations, 
the "low" fit range results are probably more trustworthy (and in any case, 
agree very well with the "normal" fit range results). 

Lastly, we have also considered the so-called "modification of the logs" 
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Prediction 


A/MeV 


d/GeV 


X 2 /dof. 


NLO ECH 
NLO MSPS 


98,116,118 
446,463,517 


1.02,0.85,0.89 
1.47,1.42,1.28 


24/46,59/88,27/44 
26/46,57/88,23/44 


LL ECH 
LL MSPS 


101,119,108 
371,417,478 


0.80,0.63,0.88 
1.12,1.02,0.83 


23/46,63/88,30/44 
24/46,49/88,22/44 


NLL ECH 

NLL ECH (mod) 

NLL MSPS 


103,107,93 
104,110,95 
200,233,268 


0.50,0.52,0.93 
0.50,0.47,0.90 
0.94,0.79,0.60 


23/46,84/88,34/44 
23/46, 81/88, 34/44 
23/46, 49/88, 23/44 



Table 1. Sensitivity of our fit values for Ajfg and C\ to the choice of fit 
range for thrust. The fit ranges are 0.05M Z /Q - 0.1 (low), 0.05M Z /Q - 0.18 
(normal), 0.11M Z /Q - 0.18 (high). 



Prediction 


A/MeV 


Ci/GeV 


X 2 /dof. 


NLO ECH 
NLO MSPS 


120,114,142 
236,115,221 


1.19,1.22,0.94 
1.65,1.84,1.27 


19/42,50/84,29/41 
19/42,68/84,37/41 


LL ECH 
LL MSPS 


124,123,128 
185,132,146 


1.06,1.07,1.01 
1.39,1.57,1.33 


21/42,51/84,29/41 
19/42,63/84,36/41 


NLL ECH 

NLL ECH (mod) 

NLL MSPS 


125,127,122 
126,128,122 
114,82,67 


0.99,0.97,1.04 
0.98,0.97,1.05 
1.29,1.49,1.62 


21/42,51/84,29/41 
21/42, 51/84, 29/41 
19/42,64/84,37/41 



Table 2. Sensitivity of our fit values for Ajj-g and G\ to the choice of fit range 
for heavy-jet mass. The fit ranges are 0.05M^/Q — 0.12 (low), 0.05Mz/Q — 
0.24 (normal), 0.125M Z /Q - 0.24 (high). 

that is often invoked in studies of event shape variables. This consists of 
modifying L = \og(l/y) — > log((2y max — y)/y) to ensure that the resummed 
parts of the expression vanish at the upper kinematic limit y ma x (which is 0.5 
for both T and pt)- The change in central values is shown in Tables 1 and 2. 
One finds that the fitted values change very little. This is to be expected since 
the restricted fit range automatically ensures that the logarithm is essentially 
unchanged in that region. 

An alternative to the approach followed here would be to apply our re- 
summations to an effective charge associated to the distribution itself similar 
to the one considered in Ref. 4J but taking the bin width to (the NLO 
approximation for this effective charge is shown as the dotted curve on Fig. 
4). This is certainly possible and leads to a p function where the LL, NLL 
and NNLL terms are known. However this creates problems with the resum- 
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Ans/MeV 




1 Ci/GeV 

0.25 0.5 0.75 1 1.25 1.5 1.75 



Figure 6: Fits to the thrust distribution for Ay^ and C\. Solid ellipses use 
ECH, dashed ellipses MSPS (with the arrows showing the effect on the central 
value of varying Q/2 < < 2Q). The ellipses indicate 2a errors generated 
by allowing \ 2 to vary within 4 of its minimum. For a summary of the data 
used see Table 3. 
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Figure 8: Examples of our best fit NLL ECH 1-thrust distributions. The solid 
curve is for Q = 44GeV, dashed is Q = 91.2GeV and dotted is Q = 183GeV. 
These predictions are compared to data from the JADE and DELPHI col- 
laborations at these energies; crosses are JADE data at 44GeV, boxes are 
DELPHI data at 91.2GeV and open circles are DELPHI data at 183GeV. 
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Experiment 


Q 


Range 


Data Points 


Source 


ALEPH 


91.2 


0.05-0.18 


7 


EU 




133 


0.04-0.15 


4 


122 


DELPHI 


91.2 


0.05-0.18 


10 


u 




133 


0.04-0.18 


5 


123] 




161 


0.04-0.18 


5 


123] 




172 


0.04-0.18 


5 


m 




183 


0.03-0.18 


11 


[23 


JADE 


35 


0.14-0.18 


2 


EH 




44 


0.12-0.18 


3 


[21 


L3 


91.2 


0.065-0.175 


4 


123 




189 


0.025-0.175 


6 


126 


OPAL 


161 


0.03-0.15 


6 


EE] 




172 


0.03-0.15 


6 






183 


0.03-0.15 


6 


m 




189 


0.03-0.15 


6 


[29 


SLD 


91.2 


0.06-0.16 


3 




30 j 


TASSO 


44 


0.12-0.16 


1 


[31 



Table 3. Summary of the data used in our fits for thrust. 



Experiment 


Q 


Range 


Data Points 


Source 


DELPHI 


91.2 


0.05-0.2 


8 


□ 




161 


0.04-0.2 


4 




23 






133 


0.04-0.2 


4 


m 




172 


0.04-0.2 


4 


[23] 




183 


0.03-0.24 


10 




23 




SLD 


91.2 


0.08-0.24 


3 




30 




ALEPH 


91.2 


0.05-0.2 


8 




21 




L3 


91.2 


0.051-0.216 


7 




2S 






189 


0.03-0.24 


14 




2E 




OPAL 


91.2 


0.0625-0.2025 


4 




21 






161 


0.0289-0.2025 


5 




2H 






172 


0.0289-0.2025 


5 




2H 






183 


0.0289-0.2025 


5 
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189 


0.0289-0.2025 


5 




n 





Table 4. Summary of the data used in our fits for heavy-jet mass. 
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mation which prevent this approach from yieding useful predictions. To see 
this, consider p LL . As explained in Sectional this is related to the leading 
log p via 

Pll(K) = p LL (TZ) + bcn 3 . (28) 

In this (unexponentiated) case is found using the one-loop beta-function 
and the double logarithmic distribution 

lda = d ex ^_ kL 2 a ^ = (2ke~ L L)a exp(-kL 2 a) = (2ke- L L)K , (29) 
a ay ay 

where k is a constant (4/3 for thrust and heavy-jet mass). This distribution 
has a peak as a function of a at a max = 1/kL 2 , and so its inverse only exists 
for 71 < 7Z(a max ) = e~ 1 /kL 2 . As a consequence Pll(JZ) vanishes at this point 
(where a branch cut starts). Adding the bcR 3 term to give pll, and later 
adding the NLL terms, does not remove this branch cut. As TZ is evolved 
from Q = oo it increases until it reaches this maximum value, and then its 
evolution becomes undefined. This value turns out to be too small to allow 
fits to the data. One could possibly "switch branches" of p at this point and 
allow TZ to decrease again, although this would of course still not provide 
a good fit to the data. Note also that this zero of p does not correspond 
to an "infrared freezing" type behaviour because p approaches the zero as 
(TZ — TZmax) 1 with 7 < 1 - thus the singularity in Eq. (jHJ) is integrable and 
the zero is reached after a finite amount of evolution in Q. 



5 Discussion and Conclusions 

In Section El we showed how it was possible to build a resummation of large 
infra-red logarithms into ECH. In principle, this allows the use of ECH for 
multi-scale observables, even when the ratio between the scales grows large. 
The approach taken was to apply a resummation at the level of the effective 
charge beta-function p. A method was described allowing such a resummed 
p to be extracted numerically to any desired accuracy from the resummed 
expression for an observable in MSPS. 

In Section|Uwe used the results of Section|3]to extend the direct extraction 
of A.j^g from e + e~ event shape observables of Ref. [1] to include a resumma- 
tion of large infra-red logarithms (in this case L = log(l/?/) where y is the 
shape variable). One could relate the observable R y (y') to an effective charge 
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1Z by exponentiation, R y (y') = exp(r TZ). Using the NLO approximation for 
1Z lead to a good fit to data. One could then numerically construct PllCR) 
and Pnll(JV) functions by resumming to all-orders the corresponding pieces 
of the p n in Eq. @. The LL and NLL predictions for R y (y') for a given 
value of Ajjg then follow on inserting these p(JZ) functions in Eq. (jUJ) and 
numerically solving the transcendental equation. To model 1/Q power cor- 
rections we fitted to a shifted distribution R PC {y) = Rpt{v — Cx/Q)- Whilst 
in principle straightforward a number of complications arose. In particular as 
1 — T approaches 1/3 the leading coefficient r goes to zero, invalidating the 
effective charge approach. This places a rather stringent upper limit on the 
fit range. There are also problems in the two-jet region due to the growing 
1Z running into a branch cut in p which appears as the image of a branch cut 
in gi in the MSPS approach. We also noted that one cannot directly relate 
the observables to an unexponentiated effective charge, as in Ref . jlj , since in 
that case Pll(R) has a different branch cut such that the energy evolution 
of 71 becomes undefined and we are unable to fit the data. Simultaneous fits 
for Ajjg and C\ were performed using data for thrust and heavy-jet mass 
distributions over a wide range of energies (see Tables 3 and 4). The 2a 
error contours in Ajjs an d C\ are shown in Figs. 6 and 7. NLO, LL and 
NLL results are shown for both standard MSPS, and for ECH. For MSPS 
there is a strong decrease in going from NLO to LL to NLL, whereas 
for ECH the fitted Ajjg values are remarkably stable. The fitted value of C\ 
is also somewhat smaller for ECH. We also investigated the stability of the 
fits to changing the fit range in Tables 1 and 2. The ECH results show more 
stability than MSPS. This, along with the stability of the fitted values with 
respect to the order of approximation, leads us to believe that the fit ranges 
we have chosen are restrictive enough to avoid the problems that appear in 
the 2-jet region for LL and NLL ECH, whilst hopefully retaining the benefit 
of including some RG-predictable sub-leading logs into our predictions. The 
fits all produce Aj^g values somewhat smaller than the world average. It is 
interesting to note, however, that they are of the same magnitude as those 
found in 4|. Converting our NLL best fit A values to a^(Mz) using the 
2-loop beta-function gives ajj^(Mz) = 0.106 for Aj^g- = lOOMeV (thrust) 
and ajfis(M z ) = 0.109 for Aj^g = 125MeV (heavy-jet mass). Similarly small 
values of a^g have also been reported in the DGE approach (see [18. ). It 
is possible that sizeable NNLO corrections (omitted by both the DGE and 
ECH resummations) might be responsible; it will be interesting to see the 
effect of matching to fixed-order NNLO results when they become available. 
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We would conclude that, notwithstanding the limited fit range and tech- 
nical complications from which ECH suffers, there is evidence that the fitted 
power corrections are reduced relative to the standard approach, although 
not as dramatically as in the DELPHI fits of Ref.jS] which are consistent 
with zero power corrections. However in that analysis corrections for bottom 
quark mass effects were made, which were not included in our analysis. Event 
shape means have also been measured in DIS at HERA [35] , and it would be 
interesting to perform an ECH analysis in that well. Unfortunately, 

in this case it is not known how to construct an ECH approximation for DIS 
where one has a convolution of parton distributions and hard scattering QCD 
cross sections. One can, however, apply PMS [S] to choose the renormali- 
sation and factorisation scales. An analysis along these lines is in progress 

EH- 
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